Single-hole properties in the t-J and strong-coupling models 
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We report numerical results for the single-hole properties in the t-J model and the strong-coupling 
approximation to the Hubbard model in two dimensions. Using the hopping basis with over 10 6 
states we discuss (for an infinite system) the bandwidth, the leading Fourier coefficients in the 
dispersion, the band masses, and the spin-spin correlations near the hole. We compare our results 
with those obtained by other methods. The band minimum is found to be at (7r/2, 7r/2) for the t-J 
model for 0.1 < t/J < 10, and for the strong-coupling model for 1 < t/J < 10. The bandwidth in 
both models is approximately 2 J at large t/J, in rough agreement with loop-expansion results but in 
disagreement with other results. The strong-coupling bandwidth for t/J > 6 can be obtained from 
the t-J model by treating the three-site terms in first-order perturbation theory. The dispersion 
along the magnetic zone face is flat, giving a large parallel/perpendicular band mass ratio. 



PACS numbers: 71.27.+a, 74.25.Jb, 75.10.Jm 
I. INTRODUCTION 

Anderson's suggestion© that the copperjpxygen planes 
of the high-temperature superconductorsH are strongly 
correlated systems has sparked renewed interest in the 
two-dimensional Hubbard model. Much of our under- 
standing of the strong-coupling limit of the model, and 
the related t-J model, has been obtained by numeri- 
cal work (reviewed in Ref. |J). Although the single- 
hole properties have been studied extensively, exact- 
diagonalization studies of small systems are hindered by 
large finite-size effects in the parameter region of interest, 
and Monte-Carlo studies of larger systems are hindered 
by the minus-sign problem; other methods have also been 
used, but there is still no general agreement on these 
properties, particularly for t/J values in the physical re- 
gion. For this reason, we have studied the sinsde-hole 
properties using the hopping basis of Trugmarotj and 
compared them with results obtained by other methods. 

In the limit U ^> t, the Hubbard Hamiltonian caa 
be approximated by the strong-coupling Hamiltoniantj 
H sc = H t -,j + H3; this differs from the t-J Hamiltonian 
H t -j (which has its own justification) by the three-site 
terms in H3: 

H t .j = -t ^2 4+s,a c i,a + J Y.( Si ' S i ~ \ n i n i) > W 
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here sites i + 8 and i + 5' are distinct nearest neighbors 
of site i, (ij) are nearest-neighbor pairs, and J = 4t 2 /U. 
The t-J and strong-coupling Hamiltonians operate in the 
reduced Hilbert space with no doubly occupied sites; this 
restriction is implicit in the above. Validity of the strong- 
coupling approximation requires U t; the parameter 
range believed appropriate to the high-temperature su- 
perconductors is 2 < t/J < 10, or 8 < U/t < 40. We 



present results for the t-J model in the region 0.1 < 
t/J < 10 and for the strong-coupling model in the re- 
gion 1 < t/J < 10. 

The single-hole properties in the t-J and strong- 
coupling models have been studied previously, the first 
having received more attention. Methods-include exact- 
diagonalization studies of small latticesOjljj-sjtiidies of in- 
finite lattices usipg-a restricted basis sculfl; Ljfc La, Monte- 
Carlo methodsE3~c3, and other methodscSEJ. Proper- 
ties discussed include the ground-state energy, the band- 
width, the dispersion, the band masses, the nearest- 
neighbor spin-spin correlations and the spectral function. 
As well, there is an extensive literature on the Hubbard 
model itself, .including recent finite-temperature Monte- 
Carlo resultsErel 

This paper studies the one-hole properties on an infi- 
nite lattice, using a restricted basis set (in effect a varia- 
tional method). Section [0] describes the basis, and Sec- 
tions IV -|V] give results for the bandwidth, the disper- 
sion, the band masses, and the nearest-neighbor spin- 
spin correlations respectively. For both models, the band 
minimum is at k = (tt/2,tv/2) and the maximum at 
k = (0,0) for the t/J values investigated. The band- 
width is approximately 2J. a . t lar ge t/J, in agreement 
with loop-expansion resultsc3E3t3 and in disagreement 
with variational Monte-Carlo resultsEil. At large t/J, the 
effects of three-site terms on the bandwidth are well de- 
scribed by first-order perturbation theory using the t-J 
ground-state wavefunction; that is, the three-site terms 
appear to have little effect on the ground-state wavefunc- 
tion at large t/J. The band mass parallel to the zone face 
is much larger than the perpendicular mass. The spin- 
spin correlations near the hole are reduced relative to the 
starting state, but remain antiferromagnetic. 



II. HOPPING BASIS 

We study a system of N — 1 electrons on a square lat- 
tice of N sites with periodic boundary conditions; the 
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Hilbert space is restricted to the S z = 1/2 sector with no 
doubly occupied sites. We use the, same basis for both 
models, nanaelv the hopping basisQcl which has been used 
previouslyD'Lnl3il3. This method allows the study of infi- 
nite systems (eliminating finite-size effects), but only cer- 
tain properties, like the bandwidth and the band masses, 
can be studied. 

In zeroth order, the basis (denoted Bq) consists of a 
single state (denoted \cN)), the Neel state with a missing 
down-spin electron. Higher-order bases are generated by 
repeatedly applying the t term in the Hamiltonian (which 
hops the hole to a nearest-neighbor site). The first-order 
basis B\ contains the \cN) state plus the four states gen- 
erated by hopping the hole. The n-th order basis B n 
consists of the states in the basis B n _\ plus those gen- 
erated by applying the hopping operator to the states in 
the difference B n ^\ — B n -2- The basis size (values are 
given in Table ||) grows exponentially with order. The 
hopping basis, which emphasizes states differing from the 
\cN) state only near the hole, cannot give a good value 
of the ground-state energy (because, for example, it does 
not generate spin interchanges far from the hole in rea- 
sonable order); the expectation is that it describes well 
properties like the dispersion and the nearest-neighbor 
spin-spin correlations near the hole. 

We have used the bases from Bq to B13 for most 
quantities, going to such large bases because some prop- 
erties were still changing significantly; even with basis 
.B13 (~ 2 x 10 6 states), however, some properties are 
incompletely converged. Various extrapolation schemes 
were considered but judged unreliable, and so we usually 
present values for the three largest bases to provide an 
estimate of the error due to the truncation of the basis. 

The system size (16 x 16; the lattice constant a is 
unity) is effectively infinite since there are no paths which 
wrap around the system even in 13-th order. Since the 
hole moves in an antiferromagnetic background, the Bril- 
louin zone is reduced to the square formed by the points 
(±7r,0) and (0,±7r). The symmetries of the lattice re- 
duce the independent part of the Brillouin zone to the 
triangle with corners at (0,0), (tt,0), and (ir/2,ir/2), de- 
noted r, M, and S respectively. Each state \n) in the 
basis is a Bloch state, an eigenstate of the translation 
operator corresponding to an allowed value of the mo- 
mentum. For each basis, and each value of the momen- 
tum k, the lowest eigenvalue and eigenvector were found 
using a conjugate-gradient method to minimize the func- 
tion (^lifl*!/) / (^\^) with respect to the expansion coeffi- 
cients in 1^) = a n \n); this method is repotted to converge 
more rapidly than others commonly usedE3, but gives the 
eigenvector to only single precision. Where necessary, the 
eigenvector was improved by a Lanczos method. 

The dispersion (in the energy as a function of JLi results 
from several processes. The Trugman pathstm trans- 
late the hole to a next-nearest-neighbor site or a third- 
nearest-neighbor site on the same sub-lattice, restoring 
the original configuration. In the lowest-order path, the 
hole hops six times around the smallest square to a next- 



nearest- neighbor site; as a result, matrix elements like 
(B2\c\ +S ulBs) are momentum-dependent. Momen- 
tum dependence can also arise from the J term in H; 
for example, the basis Bi contains states with the hole 
translated by 2a and a pair of flipped spins, and so matrix 
elements like (Bo\Si ■ Sj | i?2 ) depend on k. The results 
show odd-even effects in the order of the basis; as the 
basis size increases, Trugman paths of higher order, and 
also states differing from the starting state by nearest- 
neighbor spin interchanges, are generated. 

Related bases were also studied, in an effort to de- 
termine which states are important for the hole proper- 
ties. The hopping basis can be described symbolically as 
B n = X)fc=o h k \cN) where h is the hole hopping opera- 
tor. We define also operators S$, S12 and S20; the first 
scrambles the 8 spins at distances a and \[2a from the 
hole (giving 70 states when operating on the \cN) state), 
the second these spins plus the four at distance 2a, and 
the third the 20 spins inside a 5 x 5 square minus the 
four corner sites. If hole properties like the bandwidth 
are determined primarily by configurations which differ 
from the \cN) state only near the hole, then the bases 
£ fe /i fc S m |ciV), or (likely better) S m £ fc h k \cN), should 
converge more rapidly than the hopping basis; we find 
the opposite: when the bandwidth is plotted against 
the inverse of the log of the basis size, these modified 
bases behave like the hopping basis, except that proper- 
ties are shifted toward larger basis sizes. We considered 
also two other bases, both of which reduce the impor- 
tance of string states (in which the hole wanders without 
looping): (i) the basis ^ fc M m ft, fc |a/V) where the opera- 
tor M m removes states in which the Manhattan displace- 
ment (| a; I + \y\) of the hole relative to its initial position 
is greater than ma, and (ii) the basis J2^L (I n h) k \cN} , 
where the operator /„ removes states with more than n 
"bad bonds" (that is, it filters states according to their 
Ising energy relative to the \cN) state; the limit 00 means 
that the hop-filter combination is applied until the basis 
no longer grows, for given n). Neither the Manhattan 
nor the Ising filters improved the convergence. We con- 
clude from these numerical experiments that the single- 
hole properties are determined not so much by the spin 
configurations near the hole as by loop and string paths. 
It appears that the hopping basis, whether in its original 
form or in the modified forms we have investigated, is 
capable of only limited accuracy even if carried to very 
high order. 

III. BANDWIDTH 

Because the lattice is effectively infinite, the lowest 
energy can be found for any k. For both models, we 
found E{k) at 81 independent k values of the form 
(2im/L, 2nm/L) with n and m integers and L = 32, 
for t/J values in the range 0.1 < t/J < 10 for the t-J 
model and in the range 1 < t/J < 10 for the strong- 
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coupling model (for which the lower values of tj J are of 
little interest). 

For the t-J model, the energy is a minimum at k = S 
(and a maximum at T) for 0.1 < t/J < 10, for all bases 
used (Bq to -B13), in agreement with all previous work. 

For the strong-coupling model, the energy is also a 
minimum at k — S (and a maximum at T) for all t/J 
in the range 1.0 < t/J < 10, but only for the largest 
bases at small t/J; this result disagrees with predictions 
(based oh fits to exact-diagonalization results for small 
systemsES) that the minimum is at M for t/J < 5. For 
the smaller bases, particularly for the smaller values of 
t/J, the minimum can be at M or elsewhere along the 
zone face; for example, the minimum is at S only in 11-th 
order and higher for t/J= 1. 

Figure |l| plots the bandwidth W = E(T) - E(S) for 
both models as found using the bases Sn, -B12, and B13. 
The convergence is good for the t-J model at all t/J in- 
vestigated; it is moderately good for the strong-coupling 
model at larger t/J, but worsens at smaller t/J. The 
t-J bandwidth is approximately t for t/J < 2 and 2J for 
t/J > 2, but decreases weakly at large t/J. The strong- 
coupling bandwidth is also about 2 J (though about 20% 
larger) and also decreases as t/ J increases. The hopping- 
basis results are incompletely converged, however; the 
bandwidth is still increasing with basis size, and the trend 
is greater at larger t/J. It is possible then that the slight 
decrease which we find is due to the finite size of the 
hopping basis. 

Figure [2] compares our values for the t-J bandwidth 
with those obtained by other methods; major differences 
occur in the physical region t/J > 2. The hopnmsiasis 
results agree best with loop-expansion results ^0'E3 and 
poorly with variational Monte-Carlo resultsEJ (for up. 
known reasons); the 4x4 exact-diagonalization resultsE.il 
at large t/J are unreliable due to finite-size effects. Our 
results at large tJ-J are qualitatively consistent with the 
mean-field resulted W ss 4 J for strong coupling. 

From Figure 1, the normalized bandwidth difference 
(W sc — W t -j)/ J is almost independent of t/J for t/J > 4. 
Since (H sc — H t -j)/J = H3/J has no explicit depen- 
dence on t or J, this suggests treating the three-site 
terms as a perturbation to the t-J model. The er- 
ror in the first-order result for the bandwidth differ- 
ence AW 1 = AEt(T) - AE'i(M), where AE^k) = 
(* t _, / | J ff 3 |* t _ J )(k), is less than 2% at t/J = 10 and 
t/J = 8, but is much larger at smaller t/ J (52% at t/J = 
4). Of course the estimate for the strong-coupling band- 
width itself is much better (errors are 0.3%, 0.3%, and 
11% at t/J = 10, 8, and 4). It appears then that the 
three-site terms can be treated in first order for t/J > 6. 

Further investigation revealed that the first-order es- 
timates of the energy at S are excellent; ((H sc ) t -j — 
E sc )/W sc is 0.1%, 0.09%, 0.06% and 0.04% at t/J = 10, 
8, 4, and 1 respectively; the corresponding values at T 
are 0.4%, 0.4%, 11% and 41%. For unknown reasons, at 
intermediate t/J values the three-site terms appear to af- 
fect the r ground state strongly and the M ground state 



very weakly. 



IV. DISPERSION 



The Fourier coefficients a; m defined by 

L/2 

E(k) — ai m cos lk x cos mky 

l,m=0 



(3) 



are easily obtained by inversion from the energy as a func- 
tion of k. The symmetries of the lattice give a; m = a m i , 
and ai m = for I + m odd. The independent coefficients 
are then the 81 ai m with < I < 16, < m < I, and l+m 
even. The coefficient aoo depends strongly on the order of 
the basis, as more states important for the ground-state 
energy are generated; it affects none of our results since 
we look only at quantities (like the dispersion) which de- 
pend on energy differences. 

Of the other coefficients, an and 020 (both positive) 
are the largest, with the ratio 020/011 less than about 
0.6 for both models for the range of t/J values studied. 
The remaining coefficients are less than about O.lan in 
magnitude for both models at the t/J values studied. 
Figures || and [| plot the two leading coefficients as func- 
tions of t/J for the two models. The convergence is of 
course qualitatively the same as for the bandwidth, good 
for the t-J model at all t/J and for the strong-coupling 
model for t/J > 4, but increasingly poor for the latter 
with decreasing t/J. 

At large t/J, the values 020/ J are almost independent 
of t/J, whereas the coefficients an/ J decrease with in- 
creasing t/J. The strong-coupling coefficients are larger 
than the t-J coefficients, reflecting the enhanced mobil- 
ity due to the three-site terms. Also, at larger t/J, the 
difference {0,20/ J) sc — (fl2o/ J)t-j for the two models is al- 
most independent of t/ J, as is the difference in the values 
of an/ J, for the reason discusscd,|ia Section [II. Figures 
H and |] also plot other resultsucj for the t-J Fourier 
coefficients; the agreement is as [expected from Section 
III. Recent Monte-Carlo resultscja, available only at 
t/J = 2.5, are about 25% higher than ours. 



V. BAND MASSES 

The band masses at the band minimum, which is at S 
for both models in the region 1 < t/J < 10, are defined 
in terms of the second derivatives of E(k) with respect 
to k: 



d 2 E(k) \ 



w 



The masses were obtained by calculating E(k) at addi- 
tional points near S and using finite-difference approxi- 
mations for the derivatives. Figures I and I give results 



3 



for the masses perpendicular and parallel to the zone face 
respectively, in units of the bare mass mo — fi 2 jit. The 
parallel mass is mu<i^arflfiifj.lian_the perpendicular mass, 
as found previouslyEjEaO'EjSt^l. 

The perpendicular mass is well converged for both 
models. For the t-J model, mj_/mo is almost linear in 
t/J at large t/J, but flattens out at small t/J. For the 
strong-coupling model, m±/mo is almost proportional to 
t/J; the smaller effective mass reflects again the increased 
hole mobility relative to that in the t-J model. 

The parallel mass is much more poorly converged, es- 
pecially at smaller t/J; even at t/J = 10 (the most fa- 
vorable value), the masses change by over 5% between 
the bases B i2 and B 13 . The poor convergence results be- 
cause the energies are nearly independent of k (the mass 
is large). For large t/J, though, it appears that m||/mo 
increases only weakly with t/J for both models and that 
the two models have the same parallel mass. 

Figures || and || also give the results from Ref. de- 
rived from their dispersion results (Table 11 of Ref. |29| ) 
using the free mass m = h 2 /It, rather than the effective 
masses of their Table III. The difference is due in part 
to a genuinely different dispersion, but part of it arises 
because they used only two components in the Fourier 
expansion (the parallel mass, being large, is particularly 
sensitive to small changes in the energy). 



VI. SPIN-SPIN CORRELATIONS 

Figures [?| and || show the nearest-neighbor spin-spin 
correlation (S^ • Sj) for pairs of sites i and j near the hole, 
for the t-J model and strong-coupling model respectively. 
The momentum is k = S (the band minimum), t/J — 
2.5, and the basis is B13. In the units of H 2 = I used, the 
spin-spin correlation is -0.75 for a singlet pair of spins, 
0.25 for a triplet pair, and -0.25 for a Neel pair. The 
correlations are antifcrromagnetic, and moderately less 
than in the starting state. The "cigar" pola^u $n figures 
7 and 8 is well known from other studied 
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TABLE I. Number of states in the hopping basis versus 



order of the basis. 
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FIG. 1. Bandwidth W , in units of J, for the t-J and 
strong-coupling (sc) models as functions of t/J for the three 
largest bases used. The lines merely connect the points, here 
and in following Figures. 
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FIG. 3. The leading Fourier coefficient an, in units of J, 
for the t-J and strong-coupling (sc) models as functions of 
t/J for the three largest bases used. The dot-dash lines give 
other results for the t-J model: (a) Ref. ^ and (b) Ref. 




FIG. 2. Bandwidth W , in units of J, for the t-J model as 
functions oft/ J. The solid line gives the hopping-basis results 
(for the basis B13) and the dot-dash lines other results: (a) 
Ref. ||, (b) Ref. (c) Ref. ||, (d) Ref. ||, and (e) Ref. [jj. 
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FIG. 4. The second leading Fourier coefficient 020, in units 
of J, for the t-J and strong-coupling (sc) models as functions 
of t/J for the three largest bases used. The dot-dash lines 
give other results for the t-J model: (a) Ref. and (b) Ref. 
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FIG. 5. Band mass perpendicular to the magnetic zone face 
at the band minimum k = S, in units of the free band mass 
mo = h 2 /2t, for the t-J and strong-coupling (sc) models as 
functions of t/J. The dot-dash line gives the t-J results of 
Ref. ||. 
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FIG. 6. Band mass parallel to the magnetic zone face at 
the band minimum k = S, in units of the free band mass 
mo — h 2 /2t, for the t-J and strong-coupling (sc) models as 
functions of t/J. The dot-dash line gives the t-J results of 
Ref. || 
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-0.2898 -0.2953 -0.2687 -0.3090 -0.2937 

•- -0.2898 -•- -0.2905 -•- -0.2735 -•- -0.2937 -• 

FIG. 7. Nearest-neighbor spin-spin correlations, in units 
of h 2 , for the t-J model at the band minimum k = S for 
t/J — 2.5 using the basis B13. 

•- -0.2926 -•- -0.2583 -•- -0.2793 -•- -0.2833 -• 
-0.2926 -0.2996 -0.2732 -0.2838 -0.2833 

j- -0.2996 -j- -0.1486 -•- -0.2318 -•- -0.2838 \ 
-0.2583 -0.1486 0.0000 -0.2318 -0.2793 

j- -0.2732 -j- 0.0000 -•- 0.0000 -•- -0.2732 \ 
-0.2793 -0.2318 0.0000 -0.1486 -0.2583 

j- -0.2838 -j- -0.2318 -•- -0.1486 -•- -0.2996 \ 
-0.2833 -0.2838 -0.2732 -0.2996 -0.2926 

i- -0.2833 -i- -0.2793 -•- -0.2583 -•- -0.2926 -• 

FIG. 8. Nearest-neighbor spin-spin correlations, in units of 
Ti 2 , for the strong-coupling model at the band minimum k = S 
for t/J = 2.5 using the basis B13. 
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